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Abstract. 

The present work is our first attempt to understand the role of reconnection in the pulsar 
magnetosphere. Our discussion is based on the observationaly inferred fact that, as the pul- 
sar spins down, the region of closed corotating dipolar field lines grows with time. This im- 
plies that reconnection must take place in the magnetosphere. We argue that non-dissipative 
reconnection along the equatorial current sheet allows for the continuous channeling of pul- 
sar spindown energy into particle energy, all the way from the light cylinder to the pulsar 
wind termination shock, and we propose that this effect may account for the low cr s i wck 
values inferred by observations. We present a simple model that allows us to relate the 
magnetic diffusivity in the equatorial current sheet to an observable pulsar parameter, the 
braking index n. When n ~ 1, the global structure of the magnetosphere approaches that 
of a relativistic split monopole where the pulsar spindown energy is carried by the electro- 
magnetic field. However, for values of n > 1.5, almost all field lines close inside the pulsar 
wind termination shock, and thus most of the electromagnetic pulsar spindown energy flux 
is effectively transformed into particle energy in the equatorial current sheet. 
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1. Introduction 

In the past few years, several steps have been taken that improved our under- 
standing of the electrodynamic structure of the pulsar magnetosphere ( Mich el"! 99 II 
IContopoulos, Kazanas, Fendt 1999| (hereafter CKF); IMestel 1999I |Bogovalov 1999) 
Spitkovsky 2006 1). A canonical paradigm begins to emerge, namely that a magnetized ro- 
tating neutron star with polar magnetic field B* ~ 10 13 G, mass M„ ~ 1.4M , radius r» ~ 10 km, 
and angular velocity Q. loses rotational kinetic energy at a rate 

E kin = ^M t rlQ,D. (1) 
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through electromagnetic torques in its magnetosphere. The poloidal structure of the pulsar mag- 
netosphere remains dipolar-like (closed) from the surface of the rotating central neutron star to a 
distance 

c 

^open — — ^ 

(ric is the so called 'light cylinder'). Beyond r ope „, the poloidal structure of the magnetosphere 
becomes monopole-like (open). The total amount of open poloidal magnetic flux is of the order 
of 

%, pen = 1.23 - . (3) 

?open 

The magnetic field along the open monopole-like magnetic field lines is wound backwards in the 
azimuthal direction by the neutron star rotation, and a toroidal magnetic field is generated. The 
source of the toroidal field is a poloidal electric current whose distribution along the open field 
lines is given by 

4-7T \ ^ open J 

when < *P < ^open, and zero outside. This current distribution implies the existence of a 
singular return electric current sheet flowing along the equator and along the separatrix between 
open and closed field lines. Beyond distances greater than a few times the light cylinder distance, 
the dominant field component becomes the azimuthal one. Finally, the electromagnetic spindown 
luminosity is given by 

(5) 

Note that eqs.|3E]and|5]are valid for an aligned rotator. The case of an oblique rotator remains 
under investigation ( Bogovalov 1999 Spitkovsky 2006 , |Contopoulos 2 006 1. We have introduced 
here a cylindrical system of coordinates (r, z, <p) aligned with the central rotating neutron star. We 
have also introduced the magnetic flux *P as 



B, S Vx^-,j, (6) 

where B ; , is the poloidal component of the magnetic field. Furthermore, (...) = <9(. . .)/dt. 

In developing the above paradigm, there has been almost general agreement within the as- 
trophysics community that the global structure of the pulsar magnetosphere may be described 
by the ideal MHD formalism, except for regions of finite extent such as regions of particle 
acceleration along the magnetic field (polar gaps (e.g. |Daugherty & Hard ing 1982 1, slot gaps 
( |Harding & Muslimo v 2005) , outer gaps (e.g. |Romani & Yad igaroglu 1995 1, light cylinder dissi- 
pation regions (e.g. Me stel & Shibata 1 994 1, etc.), and current sheets which, as we will see, seem 
to be of fundamental importance to the electrodynamics of the system. The above regions where 
perfect coupling between particles and the magnetic field breaks down influence the solution of 
the ideal MHD problem through the boundary conditions. We have dealt with the presence of par- 
ticle acceleration gaps in previous papers ( |Contopoulos 2005||Contopoulos & Spitkovsky 2006) . 
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We have shown there that particle acceleration along the magnetic field 'consumes' a certain 
amount of electric potential V gap ~ 10 13 V (e.g. lHibschman & Arons 2001) from the total 'avail- 
able' electromotive potential 

V = IZZL ~ 10 13 - 17 Volt (7) 
2nr lc 

that develops on the surface of the rotating neutron star, and therefore, only the remaining 
amount is available to establish the global ideal MHD magnetosphere. In general, V gap <sc 
V, except for pulsars near their death line where V gap ~ V, in which case, the presence 
of particle acceleration gaps has important implications for the neutron star spindown (see 
|Contopoulos & Spifkovsky 2 006 1. 

In the present work, we deal with the second region where ideal MHD breaks down, namely 
the global magnetospheric current sheet. In order to understand how reconnection may work 
in the pulsar magnetosphere, we begin in § 2 with the presentation of a simple magnetostatic 
analog. In § 3, we argue that reconnection does take place, and in § 4 we show that reconnection 
extends from the light cylinder to the pulsar wind termination shock, and that it modifies the 
global structure of the pulsar magnetosphere. In § 5 we discuss the significance of reconnection 
for the acceleration of the pulsar wind to the low cr s h oc k values inferred by observations. Finally, 
in § 6 we present a summary of our conclusions. 



2. A simple magnetostatic analog 

The role of current sheets in ideal MHD is to help accomodate the requirements imposed on the 
system by the specific boundary conditions of the particular problem under investigation. In order 
to illustrate our point, let us consider a simple non-relativistic axisymmetric system consisting of 
a magnetostatic dipole surrounded by an equatorial thin disk of half-thickness 

h <k r , (8) 

with a hole inside a radius r = r open . The magnetic field structure may be obtained as a solution 
of the force-free (and current-free outside the disk) problem described by 

V x B = , (9) 

or equivalently 

— -L— + —-Q (10) 

dr 2 r dr dz 2 

with boundary conditions a) dipolar field at the center, and b) horizontal field right above and 
below the thin disk, 

B r (r<r o/ , e „;z = 0) = 0, and (11) 



B z (r > r 0P en,Z = h) = B z (r > r open ;z = -h) = 



(12) 
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These particular boundary conditions determine the magnetic field structure shown in fig. 1, 
which further implies the presence of a current sheet in the azimuthal direction all along the 
surface of the infinitely thin disk. Note that, because of the symmetry of the problem, B r (r; z) = 
-B r {r; -z), B^r; z) = -B^r; -z), and B z (r; z) = B z (r; -z). 

The reader may have already noticed that the field structure obtained in the infinite conduc- 
tivity case, fig. 1, is very similar to the poloidal magnetic field structure of the CKF solution 1 . 
Moreover, one realizes that the assumption of infinite conductivity (i.e. no field reconnection) has 
been central in obtaining the CKF solution, where we required that field lines open up beyond a 
certain radius r open , and never again reconnect. 

Now, let us assume that, once the solution shown in fig. 1 is established, we somehow in- 
troduce a small amount of resistivity (or equivalently diffusivity) in the disk. Formerly open 
magnetic field lines will now close accross the disk and will slowly begin to diffuse inward. If 
the diffusion is slow enough, we may approximate the magnetic field evolution as a sequence of 
force-free equilibria outside the resistive thin disk, coupled to the magnetic field threading the 
disk. The evolution of the latter is determined by the induction equation inside the disk, 

B = -cV xE = Vx(vxB-/7VxB). (13) 

We have made use here of the generalized Ohm's law 

v x B 77 

E + = -VxB, (14) 

c c 

where v and 77 are the matter velocity and magnetic diffusivity in the disk respectively. The z- 
component of ea.H3lvields 

B z (z = 0) = -- \-rv r B z - W —^\ , (15) 



r dr\ ' dz 
or equivalently 

dB r B r (z = h) 

*F = -2nrv r B z - 2nT]r—— ~ -2nrv r B z - 2mir (16) 

dz h 

(we have made use here of the fact that the disk is thin, and hence d/dz » d/dr). When a flow 
field is present in the disk, as is the case for the pulsar magnetosphere, the magnetic flux evolution 
will be determined by the imbalance between advection and diffusion (first and second terms in 
the r.h.s. of eq . 1 1 61 respectively) . In the present section, we will only consider the simplest case 
v = 0. A radial equatorial flow with v ~ cf will be introduced in the next section. 
One may evolve the above system of equations (eqs.llOI&U6>. as follows: 

1. Start with any initial distribution B z (r > r ope „;z = 0;f = 0) in the disk (in the example 
discussed here, we take B,(r > r ope „;z = 0;f = 0) = 0). 



1 The equatorial boundary conditions, eqs. ll 1I & I12I are the same as in the CKF problem (described by 
eq. l32l below). Therefore a solution of eg. 1 101 'kills' the first part in the l.h.s. of eg. 1321 Moreover, around 
and beyond the light cylinder, the r.h.s. balances the second part in the l.h.s. of eg. 1321 
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2. Solve eq. EH to obtain the field structure above and below the disk. This also yields the 
distribution B r (r > r open ; z = h; t — 0) on the surface of the disk. 

3. Obtain the distribution B,(r > r open \ z — 0; t — t + dt) by evolving eq.^]for one time step df. 

4. Repeat steps 1 to 3 to obtain the field evolution at all times t. 

A sequence of field equilibria is shown in fig. 2, where timescales are normalized to the charac- 
teristic time r ope „h/j]. One sees that, if we allow for an infinitesimally small amount of magnetic 
diffusivity in the disk, given enough time, the field structure will approach asymptotically the 
structure of an isolated magnetostatic dipole (in that case, B r (z = h) « B r (z — 0) = 0, and thus, 
V P = 0). We emphasize that this would not have been possible, had we insisted that the disk were 
infinitely conductive, as is the case in fig. 1 . 

Up to now, we have ignored the large radial magnetic stresses that develop inside the thin disk 
in the presence of reconnection (see fig. 2). The infinite conductivity solution shown in fig. 1 is 
stress-free, and therefore it is valid even for a non-rigid disk. However, in that case the solution 
would be unstable, and any small perturbation would result in global magnetospheric oscillations. 
Any small amount of dissipation would damp those oscillations, and the field structure would 
approach asymptotically the structure of the isolated magnetostatic dipole, which is the structure 
with the lowest possible magnetic field energy (asymptotically, B oc r~ 3 for the dipole, and oc r~ 2 
for the solution shown in fig. 1 respectively). We conclude that the sequence of field equilibria 
shown in fig. 2 will take place only if the equatorial disk is 'rigid enough' to hold the radial 
magnetic stress that develops because of the presence of reconnection. 



3. An argument for reconnection 

We will now argue that reconnection does indeed take place in the pulsar magnetosphere. We 
would like to note that, although magnetic reconnection is believed to play a very important role 
in various systems of physical and astrophysical interest (flares in the solar corona, planetary 
magnetospheres, laboratory plasmas, etc.), we know very little about the physical processes that 
control reconnection in other astrophysical systems like pulsar magnetospheres and accretion 
disks. Although it is customary to assume conditions of infinite conductivity, in the case of the 
pulsar magnetosphere, we will see that this is not justified by observations. 

One thing we are certain about pulsars is that they spin down. We also believe that the spin- 
down is due to magnetospheric torques. Equating the magnetospheric energy loss (eq.[5} to the 
observed stellar spindown (eq.Q, one obtains a relation between Cl and Q, namely that 

!_i (17) 

8 M t c r 2 open 

Furthermore, in the case of the six pulsars of table 1, Q is found to be proportional to some power 
of Q 



Q oc Q" , 



(18) 
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where n, the braking index, has been observed to have values between 1 and 3. This implies that 

n 

tH <1 (19) 



He \&o 

(and in fact, in some cases, r opm jri c <k 1). We have assumed here that n remains unchanged 
throughout the pulsar's lifetime (a simplification), and that r open = r\ c at pulsar birth when Q = 
Q„. If we assune that there is no reconnection, we imply that ^ open , and consequently r open , 
remain fixed to the value they attained at pulsar birth. This is equivalent to Q. cc £2, and therefore 
n = 1, which is not what is observed. 

We conclude that in order to account for braking index values n > 1 in our standard picture of 
electromagnetic pulsar spindown, reconnection has to take place somewhere in the pulsar mag- 
netosphere outside r open . As we will see in the next section, the natural place to look for it is the 
equatorial current sheet that develops outside r open . We may obtain an estimate of reconnection, 
by calculating the rate of decrease of open magnetic flux as the pulsar spins down, 



(«-l) Q .-.sgggg 
2 open D. ~ U year 



^ open — r. ^ open ^ ~ 10 (20) 



(see table 1). We would like to emphasize here that the fact that *¥ open + has nothing to do 
with the fact that the pulsar spins down. Even in the theoretical case where some artificial source 
of energy kept Q unchanged, as long as r open < r\ c the closed line region would grow against 
the open line region at a rate determined by eq. ^] This realization allows us to relate the mag- 
netic diffusivity in the equatorial current sheet immediately outside r ope „ to the observable pulsar 
parameters, namely that 



he 2nrri c £lB r (r, z = h) 



\B Z \ 



B r (r,z = h) 



(21) 



We have introduced here a radial equatorial flow with v r « c beyond r = r open . This accounts for 
the high Lorentz factor electron-positron outflow that is expected to flow along the open magnetic 
field lines. Inside the corrotating closed line region B r (z — 0) = (see CKF), and therefore just 
outside r ope „, B r (z - h) is also very small (figure 11 in Timokhin 2006 1. Moreover, in the ideal 
MHD case (CKF), B z (z = 0) is taken equal to zero beyond r open . In the presence of reconnection 
we expect that \B Z \ <k |Z? r | in the equatorial region beyond r open (see solution in the next section), 
therefore, we may ignore the second term in the r.h.s. of eq.|^ We thus propose a simple model 
where rj/hc depends only on the braking index n through the following approximate relation, 

f = K(n - 1) , (22) 
he 

where the parameter k will be determined in § 5. This model is consistent with our understanding 
that, the further n differs from unity, the greater the role of reconnection in the magnetosphere. 

We would like to remind the reader that "a braking index that reflects the physics of the pul- 
sar braking mechanism can only be measured for the very youngest pulsars which show either 
little glitch activity or have been observed for many years, in which case the statistical effect of 
glitches can be determined" Mobbs et al. 2 004 1. This is why we have restricted our discussion 
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to the case of the six pulsars of table 1 . However, in view of the recent measurements of O for 
many more pulsars (e.g. |Johnston & Galloway 1999] IHobbs et al. 20041 . one might be tempted 
to also consider values of n very different from the ones in table 1 . In the context of our present 
model, values of n » 1 may be interpreted as a closed line region growing through very effi- 
cient (anomalous?) reconnection ~ n times faster than the dipole canonical value 2 . On the other 
hand, values of n «: 1 may be obtained when a strong stellar wind acts to open up formerly 
closed magnetic field lines (e.g. |Harding, Contopoulos & Kazanas 19 99). In this latter picture, 
reconnection is irrelevant. One may even speculate that the wide range of n values may reflect 
a long period (100-1000 years) cyclic evolution of the system between the above two extremes 
(see |Beskin, Biryukov & Karpov 2 006 ). 



4. Continuous equatorial reconnection 

In the previous section, we have only considered reconnection in the poloidal magnetic field 
component, which is generated by a corresponding toroidal electric current. The MHD rotator, 
however, develops also a toroidal magnetic field component with a corresponding poloidal 
electric current (eq.@} such that 

B*-- 2 -!™!. (23, 
cr 

The most natural place one would expect reconnection to take place is the global magnetospheric 
current sheet where the electric current density is the greatest. We remind the reader that an 
important element of the force-free ideal MHD solution is the presence of an equatorial current 
sheet (poloidal + toroidal) that extends from the surface of the star, along the separatrix between 
open and closed field lines, and along the equator. The poloidal component of this current sheet 
constitutes the return current for the globally distributed (i.e. non-singular) electric current that 
flows along the open field lines from the stellar surface to infinity. This particular electric current 
distribution and therefore the corresponding return current sheet was obtained as an eigenvalue 
of the problem by requiring that there is no magnetic field discontinuity at the mathematical 
singularity of the light cylinder. We would like to emphasize this point again: given the particular 
equatorial boundary condition that we have chosen in CKF, namely that B z (r > r open ;z — 0) = 0, 
the existence of the above global current sheet is inevitable. 

As we will now see, the latter choice of equatorial boundary condition is not 100% valid 
when we introduce nonzero magnetic diffusivity rj and half thickness h in the equatorial poloidal 
current sheet. Because of the symmetry of the problem, 

B(r;z = 0) = fi-z, (24) 

2 Other physical mechanisms such as magnetic field decay (e.g. Urp in & van Riper 1993| may also be 
responsible for higher than normal values of n. 
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and therefore, when B z (z — 0) + a radial component of the magnetospheric poloidal electric 
field E = -rQ0 xB p /c develops along the surface of the equatorial current sheet, namely 
Qr Qr 

E r (r; z = h) = B z (r; z = h)* B 7 (r; z = 0) . (25) 

c c 

For a slowly evolving magnetosphere, the 0-component of the induction equation for the equa- 
torial current sheet (eq.H3> implies that 



dE r 
~dz 



~ , (26) 



(28) 



which further implies that E r (z = h) * E r (z = 0), or equivalently 

dB^ nB^z = h) t, 2I(V) 

ClrB : * 77— * = . (27) 

az h he r 

This may be rewritten as 

^ T) 4tt/( ( P) 
dr he Qr 

for z - and r > r ope „. Note that interior to r open , B r (r;z - 0) = 0, and therefore the field 
structure is close to dipolar with 

S-S. (29) 
or r 

/(^P) is determined similarly to CKF as an eigenvalue of the solution of the force-free magneto- 
sphere outside the current sheet (see below). 

The physical meaning of eq. is that the vertical (z) component of the electromagnetic 
(Poynting) energy flux is nonzero at the surface of the equatorial current sheet, 



— ExB 

4n 



= ^E r B^(z = h) , (30) 
z 4n 



and hence electromagnetic energy is continuously channeled to the equatorial current sheet where 
it is transformed into some different (non-electromagnetic) form. As we will see in the next 
section, it is transformed into particle kinetic energy. Another way to look at eq. I27l is that field 
lines 'corotate' with the neutron star at angular velocity Q, and therefore, particles 'trapped' 
along field lines would corotate at speeds rQ (and at distances r > r\ c these would exceed the 
speed of light) unless the equatorial region is diffusive enough for the field line to 'slip through' 
the particles at a rate given by the r.h.s. of eq.|^] This is a generalization of the argument that 
led us in CKF to propose the boundary condition that B z (r > r\ c \z — 0) = 0. 

If one knows the distribution of r\jh with distance along the equatorial current sheet, eg. 1281 
together with the special relativistic generalization of eq.|5] 

(VxB)xB+% c E = 0, (31) 

or equivalently 

r 2 \{d 2x V 1 d*F d 2y ¥\ ,2 4/0F)d/eF) 

~~ T 2 

r ic> 



dr 2 r dr + dz 2 ) + r dr c 2 dT ' ^ 
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yield the structure of the pulsar magnetosphere all the way to the wind termination shock distance 
r shock ~ ^0 9 r open . We may expect that, the presence of equatorial reconnection modifies the CKF 
solution only slightly, with magnetic flux gradually reconnecting along the equator. 

In order to obtain a solution for the structure of the global pulsar magnetosphere, we will 
make the simplification that n/h remains constant at all distances, equal to the value given by 
eg. 1221 We argued that, for small enough values of n/h, the structure of the magnetosphere will 
not differ much from that of the CKF solution, and in particular, the poloidal electric current 
distribution will not differ much from eq.[4] In that case, ea.l28lvields 
2 



"P - ^open Mn-1) ^3) 



\ f open J 



for z — and r > r open . Given the equatorial distribution of m, we obtain the global magneto- 
spheric structure shown in figs. 3 & 4 for 2 different values of n. Note that we have assumed here 
that, since reconnection is proportional to the electric current density, the perfect particle-field 
coupling required in deriving eq.[32 remains valid everywhere throughout the magnetosphere, 
except in regions with very large electric current density, as is the case for the equatorial current 
sheet. 

We are thus in a position to estimate the amount of magnetic flux that will close inside 
the pulsar wind termination shock as a function of the diffusivity parameter n/hc, and more 
specifically, as a function of the difference (n - 1). As we will see in the next section, this result 
may be related to the acceleration mechanism in the pulsar wind. 

5. Equatorial reconnection and wind acceleration 

We have shown that the small amount of reconnection implied by the observations of pulsar 
spindown requires the presence of a radial electric field component along the equatorial return 
current sheet in the pulsar magnetosphere (eq.l25>. In other words, a potential drop 

V rec = (l-^^]v (34) 

\ i open J 

develops along the equatorial current sheet, from about the light cylinder to the distance of the 
pulsar wind termination shock. Here, ^(jsHock) is the amount of magnetic flux that does not close 
accross the equatorial current sheet by the time the pulsar wind reaches the termination shock. 

According to reconnection theory (e.g. lPriest & Forbes 200 0). when the value of the recon- 
nection electric field (in our case E r as given in eq.l25> is greater than a characteristic value, the 
Dreicer field defined as 
27reVlnA 

E D = y=, , (35) 

kT e 

particles that enter the current sheet get accelerated by the electric field and do not interact during 
their acceleration process. Here, n e and T e are the particle number density and temperature in the 
current sheet respectively; In A is the Coulomb logarithm. For number densities up to a few orders 
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of magnitude greater than the local Goldreich- Julian number density ( Goldreich & Julian 1969 1, 
the radial electric field inside the pulsar magnetosphere equatorial current sheet is found to be 
orders of magnitude larger than the Dreicer field, as is the case in the geomagnetic tail (see 
|Lyons & Williams 1984| for a review). In this situation, the orbits of individual particles must be 
calculated by explicit integration of their equations of motion, and the acceleration mechanism 
is essentially a finite gyroradius effect (i.e. non-MHD). According to that scenario, energy in the 
magnetic field lines that reconnect accross the current sheet may be transformed into the kinetic 
energy of the accelerated particles without dissipation into other forms of energy l |Speiser 1987] 
refers to this process as 'inertial dissipation'). This result is quite encouraging. The equatorial 
potential drop (eq.!34> is available to accelerate some particles (electrons/positrons) up to Lorentz 
factors 

T shock = ^ f 1 - ~ f 1 - ^) 10" . (36) 

in e C y t open J \ *open J 

We have obtained here a perfect particle accelerator, all the way from the light cylinder to the ter- 
mination shock. It seems thus natural to associate the particle dominated wind that is inferred 
observationaly in the case of the Crab pulsar ( Rees & Gunn 19741 IKennel & Coroniti 19841 1 
with the particles that belong to the equatorial current sheet (mainly positrons/electrons for an 
aligned/counter-aligned rotator respectively). In a realistic pulsar, the finite angle between the 
magnetic and rotation axes, and the finite scale height h of the equatorial current sheet will lead 
to a finite opening of the equatorial particle dominated pulsar wind. Note that our picture of 
gradual poloidal + toroidal magnetic field reconnection is different from the reconnection of the 
toroidal magnetic field that has been proposed as a source of the pulsar wind acceleration in a 
'stripped wind' JCoroniti 1990llMichel 1994l|Lyubarsky & Kirk2001> . 

We have started to investigate the relativistic dissipationless acceleration of particles in the 
equatorial current sheet (Contopoulos, Gontikakis & Eftymiopoulos, in preparation). This by 
itself is a very promising new area of research. The dynamics of the equatorial current sheet 
will also be studied. In particular, in analogy to the issue of the rigidity of the conducting disk 
in our simple magnetostatic analog of § 2, we need to address the issue of the stability of the 
current sheet in the presence of the radial magnetic stresses that develop accross it due to recon- 
nection. It is interesting that the ideal MHD CKF solution implies that the current sheet is elec- 
tricaly charged (the vertical component of the electric field on the surface of the current sheet, 
E z = QrB r /c, changes sign accross it), and therefore, we expect that the radial electric field E r 
that develops along it contributes to counteract the radial magnetic stress. Note that dissipation- 
less particle acceleration due to reconnection in the equatorial current sheet was first discussed in 
|Romanova, Chulsky & Lovelace 2005] Our presentation, however, differs from theirs in that we 
tried to relate the toroidal magnetic field reconnection to the global redistribution of the poloidal 
magnetic field in the pulsar magnetosphere. This may also be seen from the fact that our equa- 
torial radial electric field scales as (r/ri c y l ~ e with e <k 1, whereas in their work it drops much 
faster with distance as (r/r/ c )~ 2 - 
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We will conclude here with a discussion of the overall energy balance in the pulsar magneto- 
sphere. In our simple picture of constant rj/h, we may estimate the fraction of the total spindown 
energy flux that remains in the from of electromagnetic field energy by the time we reach the 
pulsar wind termination shock, 



LQ¥(r shock )) ^ Up""* /ffldY = 3 l ^(r shock ) 

U^open) jj"^" /(<P)d¥ 2 I ^open 



1 



1 ^{r shock) 



3 *F 



open 



\ Vopen I 



2«(n-l) 



1 + 



1" shock 



open 



2tc(n-X) 



(37) 



(e.g. Gruzinov 2005 1. The rest of the spindown energy flux is gradually channeled to the particles 
along the equatorial sheet. It is very tempting to associate the inferred values of cr s hock (defined 
as the ratio of energy flux in the electromagnetic field to the energy flux in the particles at the 
pulsar wind termination shock), with the of energy flux estimated above remaining in the form 
of electromagnetic field energy to the energy flux transferred to the particles in the equatorial 
current sheet. It is interesting that if we fix our only free parameter to the value 



k = 0.05 , 



(38) 



we are able to produce estimates for both the Crab and Vela observationaly inferred values of 
Cshock (see table 1). 



6. Summary and conclusions 

We have presented here a new picture for the global structure of the pulsar magnetosphere, where 
field lines remain closed even beyond the light cylinder (figs. 3 & 4). This is only possible if we 
take into account the effect of magnetic field reconnection accross the equatorial magnetospheric 
current sheet. We argue that observations of braking index values n > 1 allow us to infer that 
reconnection is indeed at work in the pulsar magnetosphere, and we produce an estimate of 
the reconnection time scale (eq. |20| table 1). In fact, we propose to quantify the amount of 
effective magnetic diffusivity as being proportional to the difference (« - 1) Ceqs. l21ll22l . In our 
model, the neutron star spindown energy carried in the form of electromagnetic Poynting flux is 
gradually directed towards the equator where it is channeled to the particles along the equatorial 
sheet with no radiation losses (ea.l37>. Without considering the details of dissipationless particle 
acceleration, we argue that individual particles may reach very high Lorentz factors (eq. \36i . 
Finally, we discuss the issue of the very low (r shock values inferred at the pulsar wind termination 
shock. In our picture the regions of field and particle dominated flow are physically separated, 
and we can estimate the ratio of the respective energy fluxes. If then we assume that these regions 
somehow mix at the termination shock, we are able to estimate the values of cr shock- When we 
calibrate our diffusivity parameter (eq.l22> in order to obtain a cr shock value relevant to the Crab 
pulsar wind termination shock, all (Tshock values lie in the range 10~ 2 - 10~ 3 except for the Vela 
pulsar where cr s hock is believed to have a value close to unity. Thus, the Vela pulsar is singular, 
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because it has a much higher value of cr s hock, together with a much lower value of the braking 
index n. Our present work offers a physical picture where the above two observed parameters are 
related. It is interesting that Vela is also singular in that it shows a very strong glitch activity. This 
may contribute to a more effective loading of the magnetosphere with charges, and thus explain 
the reduced effective magnetic diffusivity in the equatorial current sheet. 

We conclude that collisionless relativistic reconnection along the equatorial current sheet 
may play a fundamental role in determining both the near and large scale structure of the pulsar 
magnetosphere, and in accounting for the acceleration of the pulsar wind, therefore, it certainly 
deserves further investigation. 
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n 

he 


T open 

(10 8 cm) 


He 

(10 s cm) 


\^open 1 

(10 3 years) 


Lfield 

u shock i . , 
^panicles 


Crab 


2.515 


0.08 


1.2 


1.6 


3 


10 2 - 10~ 3 8 x 10~ 3 


Vela 


1.4 


0.02 


0.7 


4.3 


113 


1 - 10-' 7 x 10-' 


PSR 1509-58 


2.84 


0.09 


5.8 


7.2 


3 


4 x 10~ 3 


PSRB0540-69 


2.140 


0.06 


1.2 


2.4 


6 


4 x 10~ 2 


PSR J 1846-0258 


2.65 


0.08 


8.4 


15.5 


2 


8 x 10~ 3 


PSR Jl 119-6127 


2.91 


0.10 


16.5 


19.5 


3 


4 x 10~ 3 



Table 1. Characteristic parameters for the 6 pulsars with measured values of the braking index n. 
In estimating r open we assumed an initial period of 1 msec at pulsar birth. We took here r]/hc = 
0.05(n - 1), and r^ock ~ 0.1 pc. 
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Fig. 1. The poloidal magnetic field structure of a magnetic dipole surrounded by an infinitely 
conducting thin disk with a hole inside some radius r open . The disk is denoted schematically with 
a line of x's along the equator. Length units are normalized to r open . Solid lines correspond to 
magnetic flux intervals of Q.W dipole, where *P 'dipole is defined as the amount of magnetic flux of 
a pure dipole that would cross the equator at distances greater or equal to r open . The dashed line 
denotes the separatrix *P = *P o; , e „ = 1.23 v fdi P oie between open and closed field lines. There is 
a close similarity between the present solution and the poloidal structure of the CKF magneto- 
sphere (CKF). 



16 



Ioannis Contopoulos: The role of reconnection in the pulsar magnetosphere 



t=0.1 t=0.2 




Fig. 2. A typical sequence of poloidal magnetic field configurations when we introduce magnetic 
diffusivity in the equatorial disk of fig. 1. Units same as in fig. 1. Time is normalized to the 
characteristic timescale r ope „h/T], where r\ and h <K r are the disk magnetic diffusivity and half 
thickness respectively. Note the development of large radial stresses in the interior of the rigid 
diffusive disk. In the limit t — > oo, the magnetic field configuration approaches that of a pure 
magnetostatic dipole, which is a stress-free configuration. 
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Fig. 3. The modified structure of the large scale axisymmetric pulsar magnetosphere when we 
account for magnetic diffusivity in the equatorial current sheet. Units same as in fig. 1 . We set 
rj/h = 0.08c, a value that we believe is relevant for the Crab pulsar magnetosphere (see text). One 
sees the closing of the magnetosphere accross the equatorial current sheet. Note that an amount 
of poloidal magnetic flux equal to about 0.3 x ¥ open reconnects accross the equatorial current sheet 
inside about 100r; c . 
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Fig. 4. Same as fig. 3 with rj/h = 0.02c, a value that we believe is relevant for the Vela pulsar 
magnetosphere (see text). The closing of the magnetosphere accross the equatorial current sheet 
is slower when the current sheet magnetic diffusivity is smaller. Note that the amount of poloidal 
magnetic flux that reconnects inside about 100r; c is equal to about 0.1V o ^„, i.e. smaller than in 
fig. 3 



